Parameter Dependence of Transonic Accretion onto a Black Hole 
in the Presence of Bremsstrahlung and Synchrotron Cooling 
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ABSTRACT 

The topologies are obtained for accretion flows around black holes as the parameters, 
accretion rate and specific angular momentum are varied. Paczyhski-Wiita pseudo- 
Newtonian potential is used to mimic the gravitational field of a Schwarzschild 
black hole. Both stellar mass(14M ) and supermassive(lO 8 M ) black holes are 
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flow dynamics are studied. Rankine-Hugoniot conditions are checked for possible 



considered. The effects of bremsstrahlung and synchrotron cooling processes on the 



shock formation of solution branches of the topology. Stability of the branches 
and the boundary conditions of the flow are likely to decide the uniqueness of the 
solution. 

Running Title: Black Hole Accretion with Bremsstrahlung and Synchrotron cool- 
ing 

Key Words: accretion, accretion discs - black hole physics - hydrodynamics - 
methods: numerical - shock waves 
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1. INTRODUCTION 

The presence of accretion disc was theoretically predicted even before they were 
observed, as with the case of neutron stars and black holes. Disc structures are 
expected in various astrophysical systems when it has angular momentum. Most of 
the astrophysical systems have some angular momentum associated with them. Solar 
system, galactic disc and accretion disc are some examples. The flow is assumed to 
be in form of a disc, which basically means that we have an axisymmetric system. 
As the matter spirals in, its specific angular momentum is dissipated and energy is 
lost in the form of electromagnetic radiation. Hence an accretion disc becomes a 
powerful source of luminosity. 

The study of accretion processes around a black hole using the equations of 
the general theory of relativity is a mathematically tedious task. But since the 
gravitational field is quite strong, a Newtonian description for the potential is not 
satisfactory especially close to the compact object. So as a compromise between 
the two, one can devise pseudo-Newtonian potentials which can produce the salient 
features of a black hole geometry. Also upgrading a numerical code which is written 
using Newtonian potential, can be trivially done by replacing the potential with 
an appropriate pseudo-potential. Paczyhski and Wiita (1980) introduced a pseudo- 
potential for Schwarzschild black hole which produces the positions of marginally 
stable and marginally bound orbits correctly and the efficiency approximately same 
as that of the effective potential of a Schwarzschild black hole. However it should 
be realized that while the gravitational force can be obtained as derivative of the 
pseudo-potential, one cannot relate the derivative of the effective potential to the 
gravitational force. 

More refined pseudo-potentials for Schwarzschild black hole geometry are dis- 
cussed in the literature (Nowak and Wagoner 1991; Artemova, Bjornsson and Novikov 
1996). A Pseudo-Kerr potential has been suggested by Chakrabarti and Khanna 
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(1992) which is valid only for small Kerr parameter. There are some other works in 
the literature which deal extensively with the usage of pseudo-potentials and discuss 
their relative merits (e.g. Das and Sarkar 2001; Das 2002; Das 2003; Mukhopadhyay 
2002; Mukhopadhyay 2003). In this paper, we mimic the gravitational effects of a 
Schwarzschild black hole using Paczyhski-Wiita pseudo-potential. 

Simplifying approximations have to be made to keep the problem tractable. One 
has to make a judicious choice motivated by application to a physical situation. 
Zeroth-order approximation is to assume that the disc has axisymmetry and it is in 
steady state. This results in reduction of the problem from (3+l)D to 2D. The equa- 
tions can be vertically averaged (for e.g. Shakura and Sunyaev 1973; Chakrabarti 
1989; Narayan and Yi 1994), assuming that the flow to be in hydrostatic equilibrium 
in the transverse direction which results in a hybrid 1.5D model. The assumption 
of transverse equilibrium is justified if the sound crossing time in the transverse 
direction is much less than the infall time, which in turn requires the disc to be 
thin. The results of numerical simulations of Molteni, Lanzafame and Chakrabarti 
(1994) agree with vertically averaged solution of Chakrabarti (1989). This implies 
that vertical averaging is a valid approximation. 

Thin disc of constant transverse height approximation results in a ID problem 
which can be used for testing (2+l)D numerical code (e.g. Chakrabarti and Molteni 
1993; Molteni, Sponholz and Chakrabarti 1996). The thin disc approximation is 
consistent with the requirement that the flow is highly supersonic w.r.t. angular 
velocity (e.g. Frank, King and Raine 1985; Longair 1994). In this work, we 
model the disc as thin and of constant transverse height. We also assume that the 
self-gravity of the accretion disc is negligibly small and hence the gravitational force 
experienced by an accreting element is essentially due to that of the central compact 
object. The magnetic field does not enter the dynamics of the flow as it is assumed 
to be stochastic. 
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The global solution of the accretion disc equations has been posed as a challeng- 
ing numerical boundary value problem by Narayan et al.(1997) where a relaxation 
method is used. In this paper, fourth-order Runge-Kutta method (Press et al. 1992) 
is used to solve the accretion disc equations. As Runge-Kutta method is fairly easy 
to implement, it allows to perform a rigorous and systematic study of the flow 
topologies, as the parameters are varied. Unless one finds a physical context where 
either one of these techniques is superior to the other, both are equally acceptable. 



2. THE BASIC FLOW EQUATIONS 

The equations which govern the flow are the basic conservation equations of mass, 
momentum and energy. The pseudo-Newtonian potential of Paczyhski and Wiita 
(1980) is used, which mimics the gravitational field around a Schwarzschild black 
hole to a sufficient accuracy. To describe the flow, cylindrical coordinate system, 
(r,<f>,z) is used. v(v r , v<p, v z ), p and p are the velocity, density and pressure of the 
accreting matter respectively. 

Continuity equation: 



^ + V.(pv)=0. (la) 



Euler equation: 



^ + (v.V)v + ^P + Vg = 0, (16) 
at p 

where, 

g = is 

the Paczyhski- Wiita potential, G is the gravitational constant, 
M is the mass of the black hole, r g = is the Schwarzschild radius and c is the 
velocity of light. 
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Energy equation: 



^(pe) + V.(pev)+A-r = 0, (lc) 

where, e = \v r 2 + U + - p + g + \v 2 is the specific energy, U = j^zzfi is the thermal 
energy, 7 is the adiabatic index, A is the cooling term, and Y is the heating term. 



3. SOLUTION OF TRANSONIC FLOW WITH BREMSSTRAHLUNG 
AND SYNCHROTRON COOLING 

We consider steady axisymmetric flow and use the approximation that the disc is 
thin and of constant height. Hence, velocity field of the flow is two-dimensional in 
the r<p plane. The dynamical equations of the flow are the conservation equations 
of mass, momentum and energy. The cooling effects due to bremsstrahlung and 
synchrotron radiation processes are included in the energy equation. In an electron- 
proton plasma, the expression for bremsstrahlung cooling process (Lang 1980) is 
given as, 



A brems = 1.43 x 1(T 27 N e N l T 1 l 2 Z 2 g f erg cmr z s~\ 

where, 



N;Z = 



m p + m e m p 



i.e., 



A brems = 1.43 x W-^-^T^gf, (2) 
and for synchrotron cooling process (Shapiro and Teukolsky 1983) the expression is, 
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The equipartition value for the magnetic field(S), which is obtained by equating the 
magnetic energy per unit mass to the thermal energy per unit mass is given as, 



b = J ^ y P . (3) 

V77-1 



Using the ideal gas equation, we get for the temperature T, 



,2 



p k 7 A; 



Using the above expressions for B and T we get, 



16 e 2 . e 47r 9 . 1 . 9/ a 2 um D . 9 p 
A SJ/ nc = — -( -n 2 i>{- 



3 c m e c 7(7 — 1) m e c 2 7 m p ' 



i.e., 



16 e 2 e 2 47r 1 2 jj,m p 2 1 6 2 
3 c m e c 7(7 — 1) m e c l 7 m p 

where, N e and iVj are electron and ion number densities respectively, Z is the atomic 
number, m p is the mass of proton, gf is the Gaunt factor, e is the electron charge, 
m e is the electron mass and p = 0.5 for purely hydrogen gas. 

The continuity equation: 



— (pvr) = 0. (6a) 
dr 
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The radial momentum equation: 



dv I dp .. 
v— + —f- + i> ' = 0, (66) 
dr p dr 



GM 1 2 _ GM 1 A 2 
ij ~~(r-r g ) + 2 V<t> ~~( r -r g ) + 2^' 

where, — ^jjf ^ is the Paczyhski-Wiita potential and A is the specific angular mo- 
mentum of the flow. 

The energy equation: 



dv 3 d(p/p) , A brems + A sync 

V T~ + j + ^ + — = °- 6c 

dr 2 dr pv 

The above three equations 6(a-c), are rewritten in a convenient form as, 
dv ^-^/r + jg^T^ + s^ 



dr a 2 /v — v 



(7a) 



1 dp 1 dv 1 . 
--£ + -— + - = 0, (7b) 
pdr v dr r 



and, 



dv 2a da a 2 dp ,. . . 

v— + —— + —-f + ip> = 0, (7c) 
dr 7 dr pry dr 

where, ' denotes the derivative with respect to r, v — v r: j — 1.43 x 1CT 27 in cgs 



is 



units, s = ¥r — (-^—) 2 , ^ (—^>) 2 (^ 2 -) 2 — and the polytropic relation p = Kp" 1 i 

3 c y m e c' 7(7— 1) v m e c z ' v 7 7 m p jr j jr r- r 

used to obtain the sound speed a = = y^f • 

The above set of three equations 7(a-c), are solved using fourth-order Runge- 
Kutta method (Press et al. 1992) as done in Chakrabarti (1990, see also Chakrabarti 
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1996a). The error in the fourth-order Runge-Kutta method is of the order of h 5 , 
where h is the integration step size. The solution of Narayan et al.(1997) is not 
complete as it uses only one critical point. Hence the possibility of the flow forming 
a shock is eliminated a priori. 

When the flow velocity becomes equal to the sound speed, the denominator in 
the expression for ^ (equation 7a) becomes zero. For ^ to be finite, numerator 
should also become zero. At this point, called the critical point, we use 1' Hospital's 
rule to get, 



dv ~B ± yj (B 2 - AAC) 



(8) 



dr 2A 



where, 



A = 7 + 1, 




v)+T, 



and, 



r z x 2a 7r 



The expressions for Ti, T 2 and T 3 are, 



Ti = -jia 



2 

+ 7T sa 



5 



; Ji = 39 f 



1 , Hm p 
i p 2 7A; 



) 



m. 




■jip2aa 



2a-2 
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and, 




Jipa 



2a-2 



2 4 

sa 
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where, k is the Boltzmann constant. Now the transonic solution branches are ob- 
tained by integrating from the critical point. For stellar black hole, we choose the 
mass to be 14M Q , 7 = 4/3 as distinct X-ray variability of GRS 1915+105 is inter- 
preted as due to instabilities in radiation pressure dominated disc (Greiner, Cuby 
and McCaughrean 2001). The mass is chosen as 1O 8 M for super-massive black hole, 
7 = 5/3 for gas pressure dominated disc. We consider synchrotron cooling process 
of variable efficiency The efficiency of bremsstrahlung cooling is chosen as unity. 
For a chosen accretion rate and specific angular momentum(A), we scan the r-axis 
from 1.5r 9 to 1000r 9 to find if it can be a critical point. In the case of non-dissipative 
flow as in Chakrabarti (1989), for a given specific energy there are unique values for 
inner and outer critical points. When there is energy dissipation we have some sort 
of a degree of freedom in choosing the critical points, each one corresponding to a 
different specific energy. 

3.1 The nature of critical points 

To find the critical point values, the sonic point conditions have to be used along 
with other constraints. The critical point, basically could be X-type (saddle) or 
O-type (centre) (Chakrabarti 1989). The value of the sound speed a c at the critical 
point location r c is found by using the value of the accretion rate. The following is 
the 'flow chart' (in FORTRAN 77 syntax) which is used for coding. 

(i) choose the value of the mass of the black hole and polytropic index gamma 

(ii) set the efficiency factor of bremsstrahlung cooling (zefac) and synchrotron 
cooling (synfac) 
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(iii) choose the value of accretion rate and specific angular momentum 

(iv) increment r c value from 1.5r g to 1000r 9 in a loop using r c = r c + .01 *r c * *1.5 

(v) for a chosen r c value, find a c using the expression, 
accrate = 2. * pi * r c * height * a c * p c 

= 2.*pi*r c * height *a c *(a c **2/r c — si-prime) * 1 . 5/ [zeta + sync * a c * *5) 

zeta = zefac * j * gf/(m p * *2) * (mu * m p / gamma/ k) * *alpha 

sync = synfac* 16./3.*e*e/c* ((e/m e /c) * *2) * 4. * pi/ gamma/ (gamma — 1.) * 
((l/m e /c/c) * *2) * ((mu * m p / gamma) * *2)/m p 

(vi) use the bisection routine (Press et al. 1992) to find a c satisfying the above 
expression for accrate 

The plots of a c vs r c obtained using the above routine are shown in Fig. 1 and Fig. 
2. These are representative plots showing the effect of the variation of synchrotron 
cooling efficiency factor and specific angular momentum (A) for a super-massive 
black hole of mass 1O 8 M . The value of the accretion rate is chosen as unity in 
Eddington units and the value for the adiabatic index 7 is 5/3, for a gas pressure 
dominated accretion disc. For all the cases in Fig. 1, the value of A is 1.8 in the 
units of 2GM/c. The lowermost curve corresponds to the bremsstrahlung cooling 
case without any synchrotron cooling. The synchrotron cooling factor is increased 
as 0, 10~ 7 , 10 -6 , 0.5 x 10~ 5 and 10~ 5 from the bottom to the top curve. As the 
synchrotron cooling is increased, the inner critical point ceases to exist and one would 
expect a shock free solution. There is an intermediate range of r which cannot be a 
critical point. Fig. 2 considers bremsstrahlung cooling case without any synchrotron 
cooling to show the effect of varying the specific angular momentum (A) of the flow. 
The value of A for the flow, in the units of 2GM/c, is increased as 1.6, 1.65, 1.7, 1.8, 
1.9 and 2.0 from the top to the bottom curve. 
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Fig. 1: Plot of log(a c ) vs log(r c ) for accretion rate=l in Eddington units, A=1.8 in 
the units of IGMjc and different values of efficiency of synchrotron cooling. The 
efficiency factor of synchrotron process is increased as 0, 10~ 7 , 10~ 6 , 0.5 x 10~ 5 and 
10~ 5 from bottom to top curve. The efficiency of bremsstrahlung process is unity 
for all cases. 
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log(r c ) 



Fig. 2: Plot of log(a c ) vs log(r c ) for accretion rate=l in Eddington units, efficiency 
of bremsst ranking process is unity and with no synchrotron cooling. The value of 
specific angular momentum A, in the units of 2GM/c, is increased as 1.6, 1.65, 1.7, 
1.8, 1.9 and 2.0 from top to bottom curve. 
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3.2 Flow topology 

We choose the accretion rate and specific angular momentum of the flow to be 
the free parameters. The specific angular momentum is varied in the range (1.6, 
2.5), in the units of 2GM/c . It has been suggested that flows with specific angular 
momentum close to that of marginally stable and marginally bound orbit would form 
steady accretion discs, when the self-gravity of the disc is neglected (Chakrabarti 
1996b). The location of the critical point is chosen to be between 1.5r fl and 1000r r 
In the case where only bremsstrahlung cooling process is considered, the critical 
point can be any one of the five types as shown in Fig. 3. They are named as 
X or Alp or plA or 1A or x type. X looks like the English alphabet X, Alp like 
Greek alphabet alpha which opens towards infinity, plA is reflected Alp which opens 
towards the inner boundary, 1A is plA which doesn't reach the inner boundary when 
A is high and x is X which doesn't reach the inner boundary when the accretion rate 
is high. 

Figs. 4 and 5 show the topology of the flow for bremsstrahlung cooling case and 
Fig. 6 for the case when both bremsstrahlung and synchrotron cooling are present. 
The synchrotron cooling factor is chosen as 10~ 6 . Supersonic branches are obtained 
for the outer critical points and subsonic branches for the inner critical points. The 
presence of synchrotron cooling is felt more in the subsonic branches. There is an 
intermediate range of r which corresponds to O-type critical point of Chakrabarti 
(1989). When both bremsstrahlung and synchrotron cooling are considered the 
accreting matter at large radial distance is very cold and the sound speed becomes 
small. Hence X, Alp and x type critical point ceases to exist. If the assumption 
that the efficiency of synchrotron cooling is same at all radial distances is relaxed, 
then the flow could start as a subsonic flow far away. The solution remains exactly 
the same if the mass of central black hole is changed from 1O 8 M to 14M , if the 
accretion rate and A are the same in the scaled units. But there will be a difference 
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due to the value of 7 changing from 5/3 to 4/3. 




Fig. 3: Plot of log(Mach number) vs log(Radial distance) of different components, 
X, Alp, plA, 1A and x, of the topology, for bremsstrahlung cooling case. 
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Radial distance (r) 



1000 



Fig. 4: Mach number vs radial distance of X plA - Alp X topology for accretion 
rate=1.0 and A=1.8. The outermost (supersonic) branches are of X type and it 
changes to Alp type as one moves in. There is intermediate range of r which cannot 
be a critical point and it is denoted by a -. Further moving in the critical point 
changes to plA and finally to X type again. Hence this topology is named as X plA 
- Alp X. 
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Fig. 5: Mach number vs radial distance of X 1A - Alp X x topology for accretion 
rate=5.0 and A=1.9. Note that the outermost (supersonic) branch is of x type, that 
the flow is cold enough for the value of sound speed to drop and hence become highly 
supersonic as it approaches the event horizon. 
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Fig. 6: Mach number vs radial distance topology for accretion rate=1.0 and A=1.8 
in presence of synchrotron cooling. The efficiency factor of synchrotron cooling is 
10 -6 and full bremsstrahlung cooling is used. The effect of synchrotron cooling is 
felt more in subsonic branches. This is same as Fig. 4 except for the presence of 
synchrotron cooling. 
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4. PARAMETER DEPENDENCE OF SHOCK FORMATION IN THE 
FLOW 

Mach number is defined as the ratio of the flow velocity to the velocity of sound. If 
the Mach number is less than unity the flow is called subsonic and if greater than 
unity it is supersonic. Flow is called transonic if it makes a transition from subsonic 
to supersonic flow or vice versa either continuously or discontinuously. The point 
of continuous transition is called the sonic point and the location of discontinuous 
transition is called shock. 

The fluid has to be compressible for it to be supersonic. Hence a shock is not 
possible in an incompressible fluid. Any perturbation that occurs in a real flow, 
moves with the medium and relative to the medium at the speed of sound. As 
the perturbation propagates, it steepens due to non-linearity and produces a shock 
(shock is basically a sharp discontinuity). It is possible for the flow variables to 
change discontinuously without violating the conservation laws of mass, momentum 
and energy. The strength of shock is defined as the ratio of preshock and postshock 
quantities. For a weak shock the discontinuity is small. For moderately strong 
shock, the discontinuity occurs in a very thin layer. 

Rankine-Hugoniot shock is adiabatic and there is no dissipation of energy at the 
shock but there is an increase in the entropy. The velocity jumps from supersonic 
to subsonic value and the postshock temperature is higher. Sound speed and thick- 
ness of the flow increases. In an isothermal shock, the mass and momentum are 
conserved. The energy dissipated at the shock is radiated away. The temperature 
of the preshock and postshock matter is the same as the cooling process is assumed 
to be very efficient. Sound speed and thickness of flow does not change across the 
shock. For isentropic shock, mass and momentum are conserved and specific entropy 
of the preshock and postshock matter is the same. There is loss of specific energy 
which is radiated at the shock location. 
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The foremost criterion while handling fluid dynamical problems is stability. 
Whether the flow will form a shock or not should be decided by stability. Su- 
personic flows are susceptible to shocks. The conditions for an infinitesimally thin 
and non-dissipative (Rankine-Hugoniot) shock are, 

(a) Pressure balance condition, 



Pi + P1V1 =P2 + P2V2 , (9a) 



and, 



(b) Energy flux condition, which reduces to, 

where, the suffixes 1 and 2 refer to pre-shock and post-shock quantities of the flow 
respectively. 

The Rankine-Hugoniot shock conditions are not satisfied for all values of accre- 
tion rate and specific angular momentum. The region of the parameter space that 
will form accretion or wind, with or without shocks, in the case of energy dissi- 
pationless flow, has been studied in Chakrabarti (1989). The parameter space of 
accretion rate and specific angular momentum is divided into four main regions, 
namely no shock in accretion (NSA), shock in accretion (SA), shock in wind (SW) 
and no shock in wind (NSW). 

In the case of accretion, the entropy corresponding to the inner sonic point 
is higher than that of the outer sonic point and the opposite is true for a wind 
solution. In the non-dissipative flow, as in Chakrabarti (1989), for a chosen value 
of accretion rate and specific angular momentum, there are unique values for inner 
and outer critical points. When we allow the flow to have some dissipation through 
bremsstrahlung and synchrotron radiation processes, for a chosen accretion rate 
and specific angular momentum, there is a range for outer and inner critical points, 

19 



whose specific energies are different. But the specific energy at inner critical point 
is less than that of the outer, when it forms a shock. 

When Rankine-Hugoniot shock conditions, equations 9(a-b), are satisfied the flow 
could make a transition from supersonic to subsonic branch and enter the horizon 
supersonically. In the case where there is no energy dissipation, shock forms at two 
locations inbetween the inner and the outer critical points. Due to energy dissipation 
both the shocks may approach each other or one of the shocks may disappear. Also 
there are more than one outer-inner branch combination that could form a valid 
shock solution. Inner vs outer critical points for which the shocks form are shown 
in Fig. 7(a-b) and Fig. 7(c-d). The Figs. 8 and 9 show the shock location vs outer 
and inner critical point location respectively, for valid shock solutions. In some of 
the cases, especially for A=1.8, both the inner shock and outer shock are distinctly 
visible. Due to presence of synchrotron cooling, the outer shock location is seen to 
move in. The same effect is seen in the context of white dwarf (see Chanmugam, 
Langer and Shaviv 1985, Table 1). There is a range of r cin and r cout which can 
form a shock. Among the different possible outer and inner branch combinations 
that form a shock solution, the actual flow may choose one of them decided by the 
stability and the boundary conditions. 

The code which is used for finding the parameter dependence of shock formation 
in accretion flow is based on the following 'flow chart'. 

(i) the parameters for the outer branch (accrate, lambda, a cout ,r cout , zefac, 
synfac) are read from an input file 

(ii) obtain the range of r that cannot be a critical point and decide the range for 

Tcin and T 'cout 

(iii) supersonic branch of outer critical point is obtained 
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(a) A=1.65 . 
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(b) X=1.7 
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Fig. 7(a-b): An example for r C j n vs r cout that can form a shock for (a) A=1.65 
and (b) A=1.7. The upper set of points in both (a) and (b) are for the case of 
synchrotron cooling factor 10~ 6 and lower ones are for synchrotron cooling factor 0. 
The bremsstrahlung cooling factor is unity and the accretion rate=1.0 in Eddington 
units. 
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Fig. 7(c-d): An example for r C i n vs r cout that can form a shock for (c) A=1.8 
and (d) A=1.9. The upper set of points in both (c) and (d) are for the case of 
synchrotron cooling factor 10~ 6 and lower ones are for synchrotron cooling factor 0. 
The bremsstrahlung cooling factor is unity and the accretion rate=1.0 in Eddington 
units. 
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(a) A=1.65 . 
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Fig. 8(a-b): The parameters are same as that of Fig. 7(a-b). The possible shock 
locations (r s hk) for a chosen outer critical point location (r cout ) is shown. Each shock 
location corresponds to one inner critical point location (r cin ). The left set of points 
are for the case of synchrotron cooling factor 10 -6 . 
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Fig. 8(c-d): The parameters are same as that of Fig. 7(c-d). The possible shock 
locations (r^) for chosen r cout is shown. For A = 1.8 case both inner and outer 
shock locations are distinctly visible. Both inner and outer shock merge together 
as r cout moves out. Only inner shock is possible for the case of A = 1.9. Left set of 
points are for the case of synchrotron cooling factor 10~ 6 . 
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Fig. 9(a-b): The parameters are same as that of Fig. 7(a-b). For a chosen inner 
critical point location (r cin ), different shock locations (r s hk) are possible. The right 
set of points correspond to the case of synchrotron cooling factor 10~ 6 . 
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Fig. 9(c-d): The parameters are same as that of Fig. 7(c-d). Both inner and outer 
shock locations are distinctly visible for the case of A = 1.8. For A = 1.9 only inner 
shock forms. The right set of points correspond to the case of synchrotron cooling 
factor 1(T 6 . 
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(iv) subsonic branch of inner critical point is obtained 

(v) Rankine-Hugoniot shock conditions are checked to find a shock connecting 
the supersonic branch and a subsonic branch 

(vi) for a chosen r cout , r cin is incremented in a loop and all the inner branches 
that shock with a given outer branch are found 

(vii) loop for r cout by reading next set of parameters for outer branch 
5. DISCUSSION AND CONCLUSIONS 

When solutions with shock and shock free solutions are possible for a given accretion 
rate and A, the flow might choose a shock solution as supersonic flows are susceptible 
to shocks. The inner region of the disc will be hotter when it forms a shock compared 
to a shock free disc. Hence discs with shocks may be a better candidate for explaining 
the hard X-rays and 7-rays, instead of solely relying on advection effects to increase 
the temperature (Narayan 1997). The assumption that the efficiency of synchrotron 
cooling is the same at all radial distances could be relaxed for more realistic modeling 
if matter were to start as a subsonic flow far away. It will be worthwhile to write 
the dynamical equations in curved space-time using the general theory of relativity 
before embarking into additional details of realistic modeling, as the flow velocities 
are relativistic and also the inner boundary conditions can be taken care properly. 
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